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Abstract. We present results of our preliminary analysis of acoustically active X-class flare 
fC \ of September 7, 2011. We report two acoustic sources detected via acoustic holography and 

verified by finding a ridge in time-distance diagrams. We compare the directional information 
Q ' extracted from time-distance and acoustic holography, showing a good agreement in this case. 

Vh . We report that the direction where amplitude of the wave-front is the largest lies through the 

<y} ■ strong magnetic field and sunspot, suggesting that absorption of the acoustic wave power by 

Cu ' magnetic field can be ruled out as a wave anisotropy mechanism in this case. 
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O \ 1. Introduction 

Although theoretically predicted [lj flare driven sun-quakes have only been discovered relatively 
recently [2] using SOHO/MDI data. Detection and analysis of sunquakes in the previous solar 
cycle 23 HI El El and others] has been considerably restricted by the availability of suitable 
£T) ' observations, see recent review [7] for more details. Recent launch of SDO/HMI [8j combined 

with the approaching maximum of the current Solar Cycle 24 now ensures that we can observe 
^ ■ this phenomenon in new unprecedented detail due to better data coverage as well as improved 

spatial and temporal resolution. 

Sunquakes, seen as near circular ripples, propagating outward from impulsive solar flares 
along the solar surface, appear 20-60 minutes after the flare onsets. Strong anisotropy of 
the surface wave- field has often been reported [TJ HI [9]. This has normally been attributed to 
the suppression and absorption of acoustic waves by magnetic field, which has been extensively 
studied theoretically and modelled [10]. It has also been observed that, in sunquakes, the largest 
amplitudes of the generated surface wavefront are in the direction of the apparent moving source. 

2. Data and Method 

GOES Soft X-ray flare list reports the XI. 8 class flare taking place in NOAA Active Region 11283 
on September 7, 2011, starting at 22:32 UT and peaking at 22:38 UT. The flare was accompanied 
by a CME as reported by SOHO/LASCO and CaCTUS CME event catalogues. 

We use full disk SDO/HMI intensity, dopplergram and line-of-sight magnetic field data at 
45-second cadence to produce three hour-long datacubes, which are extracted by remapping 
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Figure 1. Context images taken at 22:37:22 UT: a) HMI continuum intensity saturated (black) 
at 30% of quiet sun intensity, b) line-of-sight magnetogram, scaled between ±1 kG and c) line- 
of-sight doppler velocity on the scale between ±1.5 km s -1 . The running difference images taken 
around 22:36-22:37 UT are shown in d) for HMI magnetogram (scaled between ±200 Gauss) 
and e) continuum showing variation between —4% and 20% of quiet sun intensity. The red 
contours represent the peak 6 mHz egression power at 2 and 3 factor the quiet Sun emission 
(^ 288 m 2 s -2 ). Blue stars indicate the locations of found time-distance ridges, while orange 
arcs on images at the top row indicate the direction of the detected time-distance ridges. The 
FOV and scale are the same across a)-e). The time-distance diagrams computed for the two 
sources are presented in panels f ) for the northern Source 1 with theoretical time-distance curve 
overplotted in white and g) for the southern Source 2. In time-distance plots, distance in Mm 
is plotted along x-axis, time in minutes since 22:20:07 UT is plotted along ?/-axis. 



and derotating the region of interest using Postel projection and the Snodgrass differential 
rotation rate. The spatial resolution of the remapped data is 0.04 degrees per pixel. We apply 
acoustic holography to calculate the egression power maps from observations [HJ Q21 02] along 
with computing time-distance diagram [21 SI El El • The egression power is computed for each 
integral frequency from 3 to 10 mHz, by applying 1-mHz frequency bandwidth filters to the 
data and using Green functions built for surface monochromatic point source of corresponding 
frequency using geometrical optics approach. The choice of 1 mHz, Az/, bandwidth implies 
timing uncertainty of At = 1/Az/ = 1/1000 seconds [3J. The pupil size is set from 10 to 45 Mm. 
The basic processing follows [HI [16] . 

Time-distance diagrams are computed by selecting a source location, rewriting the observed 
surface velocity signal [v] in polar coordinates relative to the source [v{r,0,t)\ and then using 
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Figure 2. Egression power snapshot measured at 6mHz taken around 22:35:07 UT is presented 
on the left. The integration regions used in deriving statistics for Figure El are presented as white 
boxes encompassing each source. This is followed by two directional egression power snapshots 
with 6q measured from x-axis counterclockwise given in the title, and 59 = n/2. See Section 
[2] for details. All axes are in Mm, field of view is the same as in context images presented in 
Figure [TJ 



the azimuthal transformation 
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where integration is normally performed over the whole circle (i.e. a = 0, (3 = 2tt). Here, we 
use m = 0, and the integration limits are chosen to represent an arc where the signal is best 
observed. Acoustic wave-packets of a sufficiently strong amplitude relative to the noise are seen 
as a time-distance ridge that largely follows a theoretical time-distance curve. 

The holography method works by using Green's function [G + (|r — r'l, t — £'),] which prescribes 
the acoustic wave propagation from a point source, to essentially "backtrack" the observed 
surface signal [^(r, £)]. This allows us to reconstruct egression images showing the subsurface 
acoustic sources and sinks: 



r r'=b r 0=27T 
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where r' and 9 are the polar coordinates describing the vector |r — r 7 | at the solar surface. Thus, 
a, b define the dimensions of the holographic pupil. Integration over 9 normally takes place over 
the whole range [0, 27r], while here, in order to quantify and compare the helioseismic directional 
information, we use the interval of the form [Oq — 59, #o + 59, ] varying #o and taking 59 = tt/2. 
The holographic method is generally more sensitive, but vulnerable to noise, variations from 
the assumed model and is susceptible to an increased possibility of false detections, which is 
normally handled by analysing the statistical significance of the detected signal ffifTT]. 



3. Results 

The context images showing the NOAA 11283 sunspots and magnetic configuration are presented 
Figure [TJ Egression power snapshot computed at 6mHz is presented in left panel of Figure 
El It shows two clear and compact acoustic sources, northern Source 1 (SI) and southern 
Source 2 (S2). We have found directional time-distance ridges for the two sources which are 
presented at the bottom right of Figure [1] with pointing information shown as orange arcs on 
top row HMI images. No ridges were found in other directions. To verify the observed acoustic 
signatures we carry out the rms analysis of the egression power, which is presented in Figure 
[3l showing statistical significance of the measured seismic signal, particularly at 5-6mHz range, 
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Figure 3. Egression rms measured at different frequencies for Source 1 (top) and Source 2 
(bottom) as indicated in Figure [2j Time is along the x-axis, averaged egression power in m 2 s -2 is 
plotted along y-axis. Blue vertical lines represent GOES flare start and end times, red horizontal 
lines show the mean integrated rms measured away from the flare peak, with yellow horizontal 
line representing the mean plus 5a variation significance level. 

for both sources. Directional holography results are presented in two rightmost panels of Figure 
[21 showing a clear separation of the sources depending on the direction. The angles are in 
good general agreement with time-distance measurements, apparently supporting the derived 
preferred direction of surface wave-front propagation. Interestingly, as can be seen from Figure 
[U for Source 2 this direction includes the strongly magnetised plasma of the main sunspot, 
suggesting that in terms of wave-front anisotropy a mechanism different from magnetic field 
suppression of acoustic wave is at play. 

Further work is needed in understanding directional results, combined with analysis of 
RHESSI and SDO AIA observations to get better understanding of the event. 
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